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Abstract 

A variety of Yang-Baxter maps are obtained from integrable multi-field equations on quad- 
graphs. A systematic framework for investigating this connection relies on the symmetry groups 
of the equations. The method is applied to lattice equations introduced by Adler and Yamilov 
and which are related to the nonlinear superposition formulae for the Backlund transformations 
of the nonlinear Schrodinger system and specific ferromagnetic models. 

1 Introduction 

The study of set-theoretic solutions of the quantum Yang-Baxter (YB) equation was originally 
suggested by Drinfeld in pQ and has attracted the interest of many researchers. Certain examples 
of such solutions had already appeared in the relevant literature by Sklyanin [2]. Set-theoretic 
solutions of the quantum YB equation, termed as classical solutions of the quantum YB equation, 
were investigated by Weinstein and Xu [3] from a geometric point of view using the theory of 
Poisson Lie groups and the associated symplectic groupoids. The algebraic theory of the subject 
was developed by Etingof, Schedler and Soloviev More recently, the dynamical theory of the 
set-theoretic solutions of the quantum YB equation was developed by Veselov in [5] (see also the 
review article [6] ) , where the notion of the associated transfer maps was introduced as a dynamical 
analogue of the transfer matrices of the theory of solvable models in statistical mechanics. Veselov 
suggested the short term YB maps when referring to set-theoretic (or classical) solutions of the YB 
equation and which we adopt in the present work. 

Let X be a set. A map R: XxX^XxX,is called YB map if it satisfies the relation 

R (2,3) R (l,3) R (l,2) = R (X,2) R (l,3) R (2,3) ^ (1) 

regarded as an equality of composite maps from X x X x X into itself. Following the usual notation, 
is meant as the identity in the second factor of X x X x X and as R in the first and third 
factors, and analogously for BP-> 2 \ R^ 2 ' 3 ^ (see e.g. [6J). In various interesting examples of YB 
maps, such as maps arising from geometric crystals [7], the set X has the structure of an algebraic 
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variety and R is a birational isomorphism. The case of CP x CP has been recently discussed 
by Adler, Bobenko and Suris |3], in relation with the classification of the so-called quadrirational 
maps. Remarkably, the resulting list of YB maps corresponds exactly to the five different ways 
of intersection of two conies. A classification for integrable equations on quad-graphs was also 
performed by the same authors in [9], where the connection between the YB relation for maps 
and the consistency property for discrete equation on quad-graphs was addressed (see concluding 
remarks of [9]). 

More recently, the connection between YB maps and integrable equations on quad-graphs has 
been investigated in [TO]. It was shown that a systematic framework for studying such a connection 
is based on the symmetry groups of the equations. The main observation is that the variables of 
certain YB maps can be chosen as invariants of the symmetry group of the corresponding lattice 
equation. Thus, it was demonstrated that YB maps can be obtained from scalar integrable lattice 
equations admitting a one-parameter group of symmetry transformations. Moreover, this idea was 
extended to scalar integrable lattice equations admitting a multi-parameter symmetry group by 
considering the extension of the equation on a multi-dimensional lattice, as well as to multi-field 
integrable lattice equations. 

In the present work, we continue the study on the connection between YB maps and multi-field 
integrable lattice equations. In section [2] we present the explicit auto-transformations of integrable 
chains introduced by Adler and Yamilov in [11] and the interpretation of these transformations as 
discrete equations on quad-graphs. Integrability aspects of such equations are discussed in section 
[3j There we stress the link between the three-dimensional consistency property of lattice equations 
on quad-graphs and the braid type relation of the considered auto-transformations. Section 0] 
deals with symmetry methods applied to lattice equations. The existence of a symmetry group of 
transformations allows one to construct YB maps from integrable lattice equations by considering 
invariants of the transformation group as YB variables. The procedure is applied to the list of 
lattice equations of section [2] and the results are presented in section [5j In section [6J we derive YB 
maps from n-vector extensions of the nonlinear Schrodinger system. We conclude in section [7] by 
making comments on various perspectives of the subject. 

2 Explicit auto-transformations of integrable chains as 2-field quad- 
graph equations 

Adler and Yamilov introduced in [11] explicit auto-transformations for integrable chains associated 
to systems of evolution partial differential equations of the form 



u t = u xx + f(u,v,u x ) , v t 



v X x + g(u,v,v x ) . 



For example the nonlinear Schrodinger system 



Ut = u xx + 2 u 2 v 



vt = ~v xx 



-2uv 2 
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is related to the chain 



Ujx = —Uj+i - ajUj - UjV j+ i , Vj x = Vj-i + atj-iVj + Uj-iVj . 

For given k € Z, this chain admits (up to the exchange of two consecutive a's) an explicit auto- 
transformation Bfc : (uj,Vj) — » (uj,Vj), which preserves the values of the fields («.,•, Vj), j £ Z \ {k} 
and the parameters ay, j € Z \ {k — 1, &;} and changes the remaining terms of the chain by 

„ O-k—l ~ Oik ~ CUk-l — Ctk / \ 

Uk = u k + - Uk-i, v k = v k -- v k +i, (2a) 

1 - u k ^iVk + i 1 - Uk-iVk+i 

5fc = a fc _i , d fe _i = a k . (2b) 

An interpretation of the above transformations can be given in terms of deformations of simple 
paths on the multidimensional orthogonal grid. This interpretation follows naturally if one starts 
from the two dimensional grid situation studied in [12] and generalizes to higher dimension inspired 
by the multi-dimensional consistency [13 [9] and the Cauchy problem on quad-graphs [T4] . 

We consider the multidimensional orthogonal grid with set of vertices Z m , edges connecting 
neighboring vertices (ni,ri2, n m ), (nx,U2, ...,n p +l, ...,n m ) and a set of m different complex 

parameters px, P2, ■ ■ ■ ■> Pm- For each (m, ri2, n p , n m ) E Z m the edge between the vertices 
(ni,ri2, ...,n m ) and (ni,ri2, ...,n p + 1, ...,n m ) is labeled by the value /3 P , so parallel edges 

are labeled with the same /3-value. A plane perpendicular to an edge at its midpoint is called 
a characteristic plane. We consider now a simple path on the grid defined by a sequence (kj) 
of vertices with the property that each characteristic plane cuts one and only one edge of the 
path. This is an instance of a Cauchy path for quad-graph equations as was analyzed in [14J. 
At the vertex Kj we put the value of the 2-field (uj , Vj ) and we assume that the value of the 
parameter a>j of the chain coincides with the /3-value of the edge (Kj,Kj + i). Then, in the generic 
case when a k ^ OLk-x, the transformation B& can be thought as a local deformation replacing the 
(uk-x,Vk-i) ° fc ~ 1 (uk,v k ) ak (u k+ x,Vk+i) angular segment of the path by (u k -i, v k -x) " fc 
(uk,Vk) " fc ~ 1 (uk+i, Vk+x) as in figureQJa) i.e. flipping one angle of the square face to the opposite. 
Notice that if a k = ot^-Xi the transformation Bj. reduces to the identity which is compatible with 
the fact that the corresponding two consecutive edges lie on the same line and there is no angle to 
flip. 

Having this picture in mind we see that the basic building block is a square plaquette on a 
2-dimensional plane (e.g. in the 1 and 2 directions) of the grid with equations (|2b|) rewritten in 
the form 

"(1,2) = u + i «(1) > «(1,2) = v ~ i v (2) , (3) 

1 - W(1)U( 2 ) 1 - U(x)V( 2 ) 

where the subscript i inside a parenthesis indicates the displacement of the fields (u, v) in the i 
direction by a unit step on the orthogonal grid (see figure QJb) ). Following the above notation we 
list below the equations from [TT] which are relevant in the present work. We denote them by £ a , 
a = 1,2,... 5, 



£a ■ «(1,2) = E a (u,U(x),V(2),Pl,P2) , U(l,2) = E a (v, V(2) , U(l) , P2 , Pi) , (4) 
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(Uk+i,Vk+i) a k {Uk,Vk) 



«fc-i 



02 



02 



(u k ,v k ) a k -i (u k -i,v k -i) (u,v) 0i {u(i),V(i)) 

(a) (b) 

Figure 1: (a): An elementary flip representing the transformation B^. (b): The corresponding 
notation of the fields and the parameters. 



where the corresponding functions on the right-hand sides, evaluated at (x, y, z, 0%, 02), are given 
in table [T] below. 



E 1 = x + (0 1 - 2 ) 



1-yz 



F = A 1 0i0 2 (x + z)(y - x) - e 

4 X { 0i 02 02(x + z)+0 1 (y-x) 



E 2 = x + (02-0i) 



y-x 
02 + yz 



E 3 = x + (0x- 2 ) 



X + z 



y + z-0i 



E R 



(0i - 02) y z + x (72 y + 71 z) 
(01 ~ 02) x + 71 y + 72 z 



Table 1: The list of functions E a (x, y, z, 0i, 02)- The additional parameters 7$, appearing in 
function £5, depend on the parameters 0i through the relations ji 2 — 2 = 6, where 5 can be fixed 
to be or 1. In function E4, e is a free parameter. 

Remark 2.1. System £\ corresponds to equations (J3j) - System £2 is related to the chain of the 
derivative nonlinear Schrodinger system [15J deformed by its first symmetry. System £3 is related 
to the chain associated to the Kaup system [16J. The remaining lattice equations £4 and £5 are 
related to chains corresponding to the Landau-Lifshits system 

u t = u xx - — — (u 2 x + P(u)) + -P'(u) , v t = -v xx + — — (v 2 x + P(-v)) + \p'(-v) , 
for the specific cases where P(u) = e (system £4) and P(u) = 5u 2 (system £5). 



3 Braid relation and the three-dimensional consistency property 

Key properties of the transformations Mj are the relations B| = 1 , (Bj Mj + i) 3 = 1 and Mj Mi = BjBj, 
for \i — j\ > 1. The first one means that each transformation Mj is an involution. The second one, 
in view of the first, yields the following braid-type relation 

Bj_|_i Bj Bj-|-i = Kj . (5) 
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A' A 1 




Figure 2: The braid relation as the consistency property for the deformations of a path on the 
multidimensional grid. 

Let us now explain the implication of the braid relation ([5]) in terms of deformations of paths 
on the grid. By applying the composite transformations appearing on each side of equation ([5]) to 
a path ABTAEZ, as shown in figure El we get two deformations. Consider first the left hand side 
of the braid relation: Mj + i flips the r corner to K, then Mj flips the A corner to A' and finally 
Bj_|_i flips the K corner to r'. For the right hand side, Mj flips the A corner to K' then Mj + \ flips 
the r corner to i~" and finally Mj flips the K' corner to A'. Thus the original path ABTAEZ is 
transformed to the same path ABT'A'EZ through two different deformations according to figure 
El viz. 

(a) : ABFAEZ ^ ABKAEZ ABKA'EZ -*±J ABT'A'EZ , 

(b) : ABFAEZ ABFK'EZ ^ ABF'K'EZ ABT'A'EZ . 

The braid relation (|5j) guarantees that the values at the vertices r' , A' are uniquely determined 
in terms of the values at the vertices B, r, A, E i.e. independent of the way we deform the path 
in order to arrive at these vertices. This is an instance of the three-dimensional consistency for 
quad-graph equations (see [131 El [DJ, as well as [ID])- The consistency property can be checked for 
an initial value configuration as in figure El with values assigned on the vertices B, r, A, E. This 
configuration is adapted to the cubic interpretation of the YB property [IS], cf also [ID] . 

4 Invariants of symmetry groups and Yang— Baxter variables 

We start with the building block of equations on quad-graphs consisting of a system of algebraic 
relations of the form 

^(/,/(i),/(2),/(i, 2) ;/3l,^) = 0, i = l,...£, (6) 

where / = (/\ / 2 , . . . , f ) and f : 1? -> C (or CP 1 ), i = 1, 2, . . . ,£. We use upper indices to denote 
the components of the fields, since we have reserved lower indices to indicate their shifted values. 
Consider a one-parameter group of transformations G acting on the domain of the dependent 
variables of a lattice equation Q, i.e. 

G:f^&{n u n 2 ,f-e), eeC, i = l,...£, 
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where nx,n 2 € Z denote the site variables of the lattice. Let denote the lattice jet space with 
coordinates (/, fj), where by fj we mean the forward shifted values of /, indexed by all unordered 
(symmetric) multi-indices J = (j 1 , j 2 , . . . j k ), 1 < j r < 2, of order k = jfcJ. The prolongation of 
the group action of G on is 

G {k) : (f , /}) -> (^(m, n 2 , /; e), ^(m, n 2 , /; e)) , i = 1, . . . I, (7) 

where $ l (1) (ni, n 2 , /; e) = $ l (ni + 1, n 2 , /(i); e), etc. 

The infinitesimal generator of the group action of G on the space of the dependent variables is 
given by the vector field 

n , 
v = '^2Q l (n 1 ,n 2 , f)dp , where Q % (n x , n 2 , f) = -j~ $ l (ni, n 2 , /; e) 

There is a one-to-one correspondence between connected groups of transformations and their as- 
sociated infinitesimal generators, since the group action is reconstructed by the flow of the vector 
field v by exponentiation 

$ l (ni,n 2 J;e) = exp(ev)/\ i = l,...n. 

The infinitesimal generator of the action of G^ on is the associated k th order prolonged vector 
field 

n k 

vW =E E Qj(ni,n 2 ,f)d fy 
i=l #J=3=0 

The transformation G is a Lie-point symmetry of the lattice equations ([6j), if it transforms any 
solution of ([6]) to another solution of the same equations. Equivalently, G is a symmetry of equations 
([6]), if the equations are not affected by the transformation (J7J. The infinitesimal criterion for a 
connected group of transformations G to be a symmetry of equations © is 

V (2) (^(/./(l),/(2),/(l,2);/?l,/9 2 )) =0, (8) 

i = 1, 2... ,n. Equations (jSJ) should hold for all solutions of equations ([6]), and thus the latter 
equations and their consequences should be taken into account. Equations © determine the most 
general infinitesimal Lie point symmetry of the system ([6]). The resulting set of infinitesimal 
generators forms a Lie algebra q from which the corresponding Lie point symmetry group G can 
be found by exponentiating the given vector fields. 

A lattice invariant under the action of G is a function / : — > C which satisfies I(g^ ■ 
(fjj)) = I(fJ.j) for all g G G and all (/,/,) G J^. For connected groups of transformations, 
a necessary and sufficient condition for a function / to be invariant under the action of G, is the 
annihilation of I by all prolonged infinitesimal generators, i.e. 

vW(J)=0. (9) 

for all v <E Q. 
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Consider a system of lattice equations of the form © (I = 2) and suppose that it admits a 
two- parameter group of symmetry transformations G generated by two vector fields {vi, V2}. The 
invariants under the action of G on j( 2 ) are obtained from the general solution of the system of 
first order partial differential equations ([9]), by using the method of characteristics. Using these 
invariant functions P, we assign now to the edges of an elementary quadrilateral the following YB 
variables 

x l = iV,f(i)), f = i\f { 2),f), p l = ^(/( 2) , 7(1,2)), q l = ^(/ ( i,2),/(i)) • (10) 

The invariants x, y, p, q, where x = (x ,x 2 ) etc, are not independent; there exist two relations 
among them 

F(x,y,u,v) =0, i = l,2, (11) 

following from the fact that the space of G-orbits on j( 2 ) is six-dimensional. Moreover, since G is 
a symmetry group of the system of equations ([6]) , it can be written in terms of the variables (I10p , 
yielding two more relations among the invariants x, y, p, q, i.e. 

V i (x,y,p,q;(3 1 ,(3 2 ) = 0, i = 1,2. (12) 

Solving the system of equations (|12p . (jllj) for p, q in terms of x, y and assuming that the solution is 
unique, we obtain a two-parameter family of maps i?(/3i,/?2) : (x,y) 1— > (p, q). The resulting map 
R satisfies automatically the parameter dependent YB relation 

^ (2 ' 3) (/? 2 , /?s) ^ (1 ' 3) , /3s) ^ (1 ' 2) , /3 2 ) = ,R (1 ' 2) , /? 2 ) ^ (1 ' 3) , /?s) ^ (2 ' 3) (/?2, /3 3 ) • (13) 

This fact follows immediately from the interpretation of the braid relation ([5]) and the initial value 
problem discussed at the end of the previous section, in terms of YB variables (see figure [3]). More 
precisely, the LHS and RHS of the braid relation ([5]) correspond to the composition of mappings 
(a), (b) respectively, as follows 

(a) : (x, y, z) — ► {x ,y ,z) — ► [x ,y ,z) — ► {x ,y ,z ) , 

f?(2,3) p(1,3) p(l,2) 

(b) : (x, y, z) — ► (x,y,z) — ► (x , y , z ) — ► {x , y , z ) . 
The preceding considerations can be casted in the form of the following: 

Proposition 4.1. If the system of discrete equations B l = 0, i = 1,2, satisfy the three dimensional 
consistency property and admits a two-parameter group of symmetry transformations, then the map 
P2) which relates the lattice invariants $10\). satisfies the parameter dependent YB relation. 

The above framework provides a method for obtaining a YB map from a three-dimensional 
consistent lattice system of the form ^ (£ = 2). However, there are circumstances where invariance 
of a system of lattice equations under a one-parameter group of transformations suffices to construct 
an associated YB map R. This is the case for all lattice equations of table where we make the 
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#(1,3) 



^(1,2) 
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#(1,3) 



^(2,3) 



(a) 



(b) 



Figure 3: Three dimensional representation of the YB relation 

lt(2) 




Figure 4: 2-field lattice equations: the values of the fields are assigned on the vertices of a cube, 
while the YB variables on the faces. 

identification / = (u,v). In order to visualize better this situation we are going to consider a cube 
with base the elementary square and place the second field on the top face as in figure HI 
Let us examine first the system E\ which, for convenience, we rewrite below 

Pi -fo 



U(1,2)=U+ - — — — U {1) 



7 (1,2) 



1 - «(1)U(2) 



; (2) • 



(14) 



1 - U(i)U( 2 ) 

Notice that the values U(2), are not involved in the equations. Clearly, equations (|14j) are 
invariant under the one-parameter group of transformations 

G : (u, v) i-> (e £ u, e~ e v) , (15) 

with corresponding infinitesimal generator 

vi = ud u - vd v . (16) 

We consider now the space J^(%)), i.e. the space obtained from J^ 2 ) by factoring out vny On 
J( 2 )(uqO we define the following lattice invariants along the orbits of vi 



(x 1 ^ 2 ) = (U/U(i) , l)U(i)) , ijp l ,V 2 ) = {U(2)/U(l )2 ) , U2«(l,2)) , 

(y 1 ,^ 2 ) = (U(2)/U, UV {2 )) , (<?W) = (tt(l,2)/W(l) , U(l,2)) 



(17a) 
(17b) 
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The space of G-orbits on j( 2 )(?;rn) is six dimensional and thus there exist exactly two functionally 
independent relations among the invariants (|17|) . namely 

pV = *V, p 1 v 2 = y 1 y 2 - (18) 

Moreover, since G is a symmetry, equations (|14p can be written in terms of the invariants (|17p . in 
the following form 

p 2 = y 2 p, q 2 = x 2 + y 2 (l-P), where P = 1 + ^ ~ ^ . (19) 

a; 1 — y z 

Solving the system of algebraic equations (|18p . (|19|) for (p,q) in terms of (x,y), we obtain the 
unique solution 

(p\p 2 ) = (V 1 P- 1 ,y 2 P), (q\q 2 ) = (x 1 P,x 2 + y 2 (l - P) ) , (20) 

which defines the YB map R : (x l ,y l ) i— > {p l ,q l ) (shaded triangles of figure S]). 

On the other hand, one could equally well consider on the space J^ 2 ) («(2) ) the following invariants 
under the action of G (blank triangles of figure H|) 

(x 3 , x 2 ) = (v/v {1) , v U( 1} ) , (p 3 ,/) = (u( 2 )/f(i,2) , U( 2 ) «(i )2 )) , (21a) 

(y 3 ,2/ 2 ) = {V(2)/v , UV {2) ) , (<? 3 ,? 2 ) = («(1,2)M, «(1)«(1,2))- ( 21b ) 

The above invariants are related by 

p 3 q 3 = x 3 y 3 , p 3 q 2 = y 3 x 2 , (22) 
while equations (|14p can be written in terms of them as follows 

/ = y 2 +x 2 (1 _ p }) (? 2 =x 2p ) where p =1 _ (fo- fo)y 3 (23) 

1 — x z y 6 

Solving the above equations for (p l , q l ), we obtain the YB map 

(p 3 ,p 2 ) = (y 3 P~ 1 ,y 2 + x 2 (l- ?)), (q 3 ,q 2 ) = (x 3 P,x 2 P). (24) 

The variables (|17|) (here called Yang-Baxter variables) were introduced in a different context in 
|llj . where the reduction of the lattice equation ([3]) to the map ()20p was also presented, together 
with an associated Lax pair. It seems that one of the motivations in [11] for this reduction was the 
illustration of the discrete Liouville theorem for integrable maps. 

We close this section by presenting a Lax pair, different from the one given in [TT], for the YB 
map (|20|) . The YB map (I20p admits the discrete zero curvature representation 

W(y\y 2 ,p 2 ,X)W(x\x 2 ,p 1 ,X) = W(p\p 2 , fa, A) W(q\ q 2 , fa, A) , 

where 

wteM=[ e i A : ? ;f {2 ). 
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A hint for considering the matrix W of the above form, is based on the YB map itself. Indeed, in 
view of equations (f2"U|) giving p 1 , one notices that 



y ^-y 2 + /3i-/? 2 ' 

can be written as a linear fractional transformation induced from the linear transformation with 
matrix W{y 1 , y 2 , Pi) on (sc , 1) (column vector), cf [T8] . 



5 YB maps arising from equations £ a 

In this section we present the YB maps obtained from the lattice equations £ a by considering 
invariants of their symmetry groups. For equations admitting a two-parameter symmetry group 
we consider as YB variables a complete set of joint invariants on j( 2 ). In the case of a one- 
parameter symmetry group, as we discussed in the previous section, one can consider two cases 
for YB variables: (i) lattice invariants on the space jt^iim) (shaded triangles in figure UJ), or (ii) 
lattice invariants on the space (blank triangles in figured]). We present YB maps obtained 

from case (i) only. For each system the analysis of case (ii) yields identical map as in case (i), up 
to appropriate rearrangement of variables. 

5.1 System £2 

Equations £2 admit the symmetry generator (|16p . thus we take as YB variables the ones defined 
by equations (|17|) . leading to the same functional relations (|18p . On the other hand, system £2 can 
be written in terms of the invariants (|17|) as 

p 2 = y 2 P, q 2 = x 2 Q, (25) 

where 

Solving the system of equations (|18|) . (|25|) for {p l ,q l ) we obtain the YB map 

(p 1 , p 2 ) = (y 1 P- 1 , y 2 P) , (q\q 2 ) = (x 1 P,x 2 Q). 

5.2 System £ 3 

Equations £3 admit the symmetry generator 

v 2 = d u - d v . (26) 

On the space j( 2 )(t>(i)) we define as YB variables the following invariants 

(x 1 , x 2 ) = (u- n (1) , v + u (1 )) , (p\p 2 ) = (u( 2 ) - U(i,2) , vq) + U(i j2 )) , (27a) 
(y\ y 2 ) = («(2) u , it + U( 2 )) , (g 1 , g 2 ) = (u(i j2 ) - , tt(i) + U(i, 2 )) • (27b) 
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They are related by 

pi + g 1 =X 1 +y\ pl +p 2=yl + y2, (28) 

and system £3 is written in terms of the invariants (I27p as follows 

p 2 = y 2 P, q 2 = x 2 P' 1 , (29) 

where 

_ y 2 - x 1 - fa 



y 2 -x 1 - Pi' 

Solving the system of equations (|28|) . (|29p for ) we obtain the YB map 

(p V) = {y 1 + y 2 (i-P), y 2 P) , (?\ q 2 ) = (x 1 -y 2 (i-P), x 2 p- 1 ) . 
5.3 System £ 4 

In the generic case e 7^ 0, the YB variables are the invariants of the symmetry group generated by 
the vector field (j26[) . System £4 is written in terms of the invariants (|27p as follows 

p 2 = y 2 -P, q 2 = x 2 + Q, (30) 

where 

1 1 (3 1 (3 2 y 2 x 1 + e 1 1 p 1 p 2 x 2 (y 2 - x 1 - x 2 ) - e 

Vi /V (hy 2 -0ixl ' W V2 /V /3 1 x2 + /? 2 ( 2/ 2_ x i_ x 2 ) • 

Inserting equations (130)) into ([28]) and solving for (p 1 ,^ 1 ) in terms of (x*,y*), we arrive at the YB 
map 

(p\p 2 ) = (y 1 + P,y 2 -P), (q\q 2 ) = (x 1 - P , x 2 + Q) . 

We examine now the case e = 0. In this case there exists, in addition, the scaling symmetry 
generated by the vector field 

V3 = ud u + vd v . (31) 
On j( 2 )(«n)) we define as YB variables the following invariants under V3 

(x\x 2 ) = (lt/«(i) , V/ U{1) ) , {p l ,p 2 ) = («(2)/«(l,2) ) «(2)/«(l,2)) . ( 32a ) 

(y 1 , y 2 ) = («(2)/« > u (2)/«) , (?\ q 2 ) = («(i,2)/«(i) > «(i,2)/«(i)) . ( 32b ) 

which are related by 

p 1 g 1 =a; 1 2/ 1 , y 2 p 1 = y 1 p 2 . (33) 
System £"4 (e = 0) can be written in terms of the invariants (|32j) . as follows 

p l =x 1 P, q 2 = x 2 Q, (34) 

where 

(/3 2 -/3i)(l + y 2 )(l-^) 



P= 1 + 
Q = 1 + 



/3 2 xi(l + y2) +/3l (i_ x i) ' 

1 (/3i-/? 2 )(l+x 2 )(xV-x 2 ) 
x 2 (3i (1 + x 2 ) + 02 {x l y 2 - x 2 ) ' 
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From relations (|33p . (|34|) we obtain the following YB map 

(p\p 2 ) = (y 1 P- 1 , y 2 P- 1 ) , (q\q 2 ) = (x 1 P , x 2 Q) . 

Finally, considering the group generated by {v2,V3}, given by (|26|) . (j3"Tj) . we take on j( 2 ) the 
following joint invariants 

(.V) = (^,^), ^.(M^,!^,, (35 a) 

y V + U {1) V + U {1) ' V «(2)+«(l,2) «(2)+«(l,2) 

feV ) = (^,^), ^.(iieazM.ieazM). (35 b) 

V M + V(2) n + V( 2 ) 7 W(l)-«(1,2) «(1)+V(1,2) 



The above invariants are functionally related by 

(l + x 1 )(l + y 1 ) (l-x 2 )(l-y 2 ) l + x 1 y 2 



(36) 



(1 + p 1 ) (1 + q 1 ) (1 - p 2 ) (1 - q 2 ) 1 + p 2 q 1 ' 
In terms of the invariants (|35p . the system £4 (e = 0) can be written as follows 

q 1 =X 1 P, p 2 =y 2 p-\ (37) 

where 

/? 1 (l-y 2 ) + /3 2j / 2 (l + x 1 ) 
ftx 1 (y 2 -l)+/3 2 (l + x 1 ) ' 

Using relations (j36j) . (j3"Tj) we obtain the YB map 



! 2 ^ (l +g l)(l+yl) \ 



(38a) 



( ^^ = ^P,(l_^Hl_^ + 1 j. (38b) 
5.4 System £ 5 

We examine first the generic case 5 = 1, i.e. 7? — /3 2 = 1. System £5 admits the scaling symmetry 
with generator V3, thus we use the invariants (|32|) . In terms of these variables the system £5 is 
written as follows 

p 1 = x 1 P , q 2 = X 2 Q, (39) 

where 

p _ (f3i-p2)y 2 + l2 + lix l y 2 Q _ (/3 2 -/?i)4^+72+7i^S 2 
(Pi -/? 2 )ic 1 +71 +~/2X 1 y 2 ' (P2 - Pi)x 2 + 71 +72X 1 y 2 

The YB map obtained from relations (|33p . (I39p is 

(p\p 2 ) = (y 1 P- 1 , y 2 P- 1 ) , (q\q 2 ) = (x 1 P , x 2 Q) . 

When 5 = 0, the system £5 becomes system £4 with e = directly (case 7, = — or indirectly 
under the gauge transformation (u,v) 1— > (— l) ni+n2 (u, u ) (case 73 = /%). 
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6 Discrete vector Schrodinger system 



In this section we show how the preceding considerations can be generalized to vector extensions 
of integrable lattice equations yielding n-vector (two-field) YB maps. This is illustrated for the 
lattice n-vector Schrodinger system 

V ' 1 - «(1) • «(2) 1 - «(1) • »(2) 

introduced by Adler }19j . where u = (u , . . . ,u n ) etc and • stands for the usual Euclidean inner 
product. Equations (|4Up represent the superposition formula of the Backlund transformation for 
the n-vector nonlinear Schrodinger system found by Manakov |20| . 

System (|40p remains invariant under the group of transformations generated by the vector fields 

Vj = u l d u i — v l d v i , 

i = 1, ... n, where no summation on indices is assumed. In order to make the presentation more 
concise, we represent the vectors u, v and their shifts in terms of diagonal matrices denoted by the 
corresponding capital letter, e.g. u — > U = diag (u 1 , u 2 , . . . , u n ). On the lattice jet space J^^i^)) 
we consider the following YB variables 

(Y\Y 2 ) = {UU {1 f\VU (1) ), (P\P 2 ) = {U (2) U (1 , 2) -\V (2) U {1 , 2) ), 
{Y\Y 2 ) = {U {2) U-\UV {2) ), {Q\Q 2 ) = {U m U {l) -\U [x) V {l , 2) ). 

The above invariants are functionally related by 

p 1 Q 1 = X 1 Y 1 , P 1 P 2 = Y 1 Y 2 , (41) 
and the lattice equations (|40|) are written in terms of them as follows 

P 2 = Y 2 S, P 2 + Q 2 = X 2 + Y 2 , (42) 

where 

S = I+ (ft - tk) (x'r 1 



l-tr(Y 2 (Y 1 )- 1 ) ' 

Solving equations (|41|) . (|42p for (P l ,Q l ) in terms of (X l ,Y l ) we obtain the YB map 

(P\P 2 ) = (Y 1 S~ 1 , Y 2 S) , (Q\Q 2 ) = (X 1 S , X 2 + Y 2 (I - S)) . (43) 

It would be interesting to investigate the relation of the YB maps studied recently in [211 [22| [23] 
with the YB map (|43l) after imposing appropriate reality conditions on the latter, as it was done 
for system £\ in [TTj . 



13 



7 Conclusions and perspectives 



We have applied symmetry methods to multi-field integrable lattice equations in order to construct 
YB maps. The main idea is that the invariants of the symmetry groups of the lattice equations 
serve as variables of the associated YB maps. 

According to [8], a birational map R : (x,y) i— > (p(x , y) , q(x , y)) is called quadrirational if 
the maps s i— > p(s, y) and t t— > q(x, t) for generic values of y, x, are invertible. This property is 
equivalent to the nondegeneracy property which is often imposed additionally in the study of YB 
maps, cf [I]. We have to notice that all YB maps we presented in this work do not satisfy the 
quadrirationality property. This can be seen immediately since {p l -,p 2 ) is independent of x 2 . Two 
additional properties which characterize certain YB maps are invertibility and unitarity. A YB 
map R(Pi, fa) with parameters {fa, fa) is called unitary if it satisfies the relation 

R {2 ^{fa,fa)R{fa,fa)=id, 

where i?( 2,1 ) = PRP and P the permutation map i.e. P : (x,y) >— > (y,x). The YB maps we 
obtained here are invertible and unitary. 

Recently, Reshetikhin and Veselov [23] studied the Hamiltonian theory of YB maps through 
Lie-Poisson structures. The Poisson structure corresponding to the auto-transformation ([3]) was 
presented in [11]. It would be interesting to explore the Poisson structure of the remaining auto- 
transformations and the associated YB maps using the theory of Poisson Lie groups as in [U [21] . 

The protagonistic role of the YB relation in quantum integrability is well-known. In view of 
recent developments (see [6] and references therein) it seems that the YB relation for maps acquires 
a prominent role in (classical) discrete integrability as well. Here, we saw that integrable lattice 
equations admitting a symmetry group can be reduced to YB maps. It remains an open question 
how integrable lattice equations which do not admit an appropriate symmetry group can be casted 
in the form of a YB map. 

By the way of studying the multi-field lattice equations introduced by Adler and Yamilov the 
question of classification of two-field lattice equations by imposing the three-dimensional consis- 
tency property is naturally raised. Before tackling this problem in full generality one can consider 
the restricted class of equations of type ([3]). It would be interesting to investigate whether the 
classification is exhausted essentially by the lattice equations in [llj. 
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